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Abstract. We present ISO observations of several H2 
pure-rotational lines (from S(0) to S(5)) towards a sam- 
ple of 16 molecular clouds distributed along the central 
~ 500 pc of the Galaxy. We also present C 18 and 13 CO 
J = 1 — > and J = 2 — > 1 observations of these sources 
made with the IRAM-30m telescope. With the CO data we 
derive H2 densities of ifj 3 - 5 ~ 40 cm~ 3 and H2 column den- 
sities of a few 10 22 cm -2 . We have corrected the H2 data 
for ~ 30 magnitudes of visual extinction using a self- 
consistent method. In every source, we find that the H 2 
emission exhibits a large temperature gradient. The S(0) 
and S(l) lines trace temperatures (T) of ~ 150 K while the 
S(4) and S(5) lines indicate temperatures of ~ 600 K. The 
warm H2 column density is typically ~ 1 — 2 x 10 22 cm~ 2 , 
and is predominantly gas with T=150 K. This is the first 
direct estimate of the total column density of the warm 
molecular gas in the Galactic center region. These warm 
H2 column densities represent a fraction of ~ 30% of the 
gas traced by the CO isotopes emission. The cooling by 
H 2 in the warm component is comparable to that by CO. 
Comparing our H 2 and CO data with available ammonia 
(NH3) observations from literature one obtains relatively 
high NH3 abundances of a few 10 -7 in both the warm and 
the cold gas. A single shock or Photo-Dissociation Region 
(PDR) cannot explain all the observed H2 lines. Alterna- 
tives for the heating mechanisms are discussed. 

Key words: ISM: clouds - ISM: molecules - ISM: dust, 
extinction - Galaxy: center - Infrared: ISM: continuum - 
Infrared: ISM: lines and bands 



region there are molecular clouds and huge cloud com- 
plexes like Sgr B2 which can be as large as 70 pc, with 
masses of 10 6 solar masses. The physical conditions in the 
GC clouds differ appreciably to those of the clouds of the 
disk of the Galaxy. The GC clouds have average densi- 
ties of ~ 10 4 cm~ 3 instead of 10 2 cm -3 typical of the 
clouds of the disk. In addition, with widespread high tem- 
peratures (up to 200 K), GC clouds are hotter than disk 
clouds. 

The temperatures of the warm gas are known mainly 
by observations of ammonia (NH3) mctastablc lines. 
Glisten et al. (1981, 1985) derived rotational temperatures 
(T ro t) of 60-120 K in several GC clouds, most of them 
in the Sgr A complex. Morris et al. ( [1983 ) showed that 
Tot ^ 30 — 60 K are common in the region \l\ < 2°. The 
most complete study of the temperature structure of the 
molecu lar ga s in the GC, was carried out by Hiittemeister 
et al. ( 1993 ). They presented a multilevel study of NH 3 
metastable lines of 36 molecular clouds distributed all 
along the "Central Molecu lar Zone" (CMZ, in notation 
of Morris & Serabyn |1996| ) and the "Clump 2" complex, 
which, although not belonging to the actual CMZ, exhibits 
similar properties. They detected warm gas at all galac- 
tic longitudes and showed that the NH 3 emission can be 
characterized by two temperature components since the 
Tot derived from the (1,1) and (2,2) levels is - 20 - 30 K 
and that derived from the (4,4) and (5,5) levels is ~ 70- 
200 K. Unfortunately, the a priori unknown abundance 
of the NH3 molecule has made it difficult to estimate the 
total column density of warm gas in the GC clouds. 



1. Introduction 

The interstellar matter in the inner few hundreds parsecs 
of the Galaxy (hereafter GC) is mainly molecular. In this 

Send offprint requests to: N.J. Rodriguez-Fernandez 

* Based on observations with ISO, an ESA project with in- 
struments funded by ESA Member States (especially the PI 
countries: France, Germany, the Netherlands and the United 
Kingdom) and with the participation of ISAS and NASA. 



The heating of the molecular gas over large regions 
(~ 10 pc) where the dust temperature is lower than 
30 K (Odenwald & Fazio gggl Cox & Laureijs [1989, 



Martin-Pintado et al. 1999, Rodriguez-Fernandez et al. 
pooo| ) is a puzzle. Indirect arguments such as the large 
widths of molecular lines or large abundances in gas phase 



of molecules such asSiO (Martin-Pintado et al. 1997, 
Hiittemeister et al. 199§|) or NH3 points towards a mechan- 
ical heating. Wilson et al. ( 1982j ) proposed the dissipation 
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of turbulence induced by differential Galactic rotation as 
a possible heating source. 

For the first time, we have measured the total col- 
umn densities of warm gas in the GC clouds by observing 
the lowest H2 pure-rotational transitions with the Infrared 



Space Observatory (ISO; Kessler et al. 1996). The H 2 pure- 



rotational lines trace gas with t emper atures of a few hun- 
dred K (see Shull & Beckwith |1982| for a review on the 
properties and the notation of the H2 molecule). ISO has 
detected H 2 pure-rotational lines in a variety of s ource s 
such as: Young Stellar Objects (V an de n Ancker 1999| ); 
galactic nuclei (see e.g. Kunze et al. 1999[ ); Photo Di ssoci- 
ation Regions (PDRs) like NGC 702 3 (Fue nte et al. |l999 , 
2000| ) or S140 (Timmermann et al. |l996| ); sh ock-ex cited 
sources such as Orion Peak 1 (Rosenthal et al. 2000); and 
proposed x-ray excited regions (XDRs) like RCW 103 (see 
Wright EOOOD . 

Our sample consists of 18 molecular clouds from 
the samples of Huttemeister et al. (1993) and Martin- 
Pintado et al. (1997). Two of these show a non-equilibrium 
H2 ortho-to-para ratio and have been studied in detail 
by Rodriguez- Fernandez et al. (2000). In this paper we 
present the other 16 clouds of the sample. The clouds are 
distributed along the CMZ, from the Sgr E region to the 
vicinity of Sgr D and the "Clump 2" complex. Four clouds 
are located in the Sgr C complex, three in the vicinity of 
Sgr A (two are in the radio Arc). Two clouds are sit uated 
in the cold dust ridge reported by Lis & Carlstrom ( 1994 ) 
that seems to connect the radio Arc and Sgr B. Other 
three clouds belong to the Sgr B complex. 

This paper is organized as follows. In Sect. 2 we present 
C 18 and 13 CO IRAM-30m observations and H 2 ISO ob- 
servations. The analysis of the CO isotopes and H2, is pre- 
sented in Sect. 3 and 4, respectively. The results and the 
possible heating mechanism of the warm gas are discussed 
in Sect. 5. 



2. Observations 

2.1. IRAM 30-m observations and results 

We have observed the J = 1 — > and J = 2 — > 1 lines of 
13 CO and C 18 with the IRAM 30-m telescope (Pico de 
Veleta, Spain) towards the GC molecular clouds given in 
Table [p. This table also gives the pointing positions and 
the complexes where the clouds belong. Figure [j] shows 
the position of the sources overlayed on the large scale 
C 18 O(l-^0) map of Dahmen et al. ( |1997| ). The observa- 
tions were carried out in May 1997, May 1998 and June 
2000. The J=1^0andJ = 2^1 lines were observed 
simultaneously, with two 512 x 1 MHz channel filter banks 
connected to two SIS receivers at 3 and 1.3 mm. The re- 
ceivers were tuned to single side band mode. The image re- 
jection, checked against standard calibration sources, was 
always larger than 10 dB. Typical system temperatures 
were ~ 250 K for the J = 1 — > lines and ~ 400 K for 



the J = 2 — > 1 lines. The velocity resolution obtained with 
this configuration was 2.7 and 1.4 km s _1 at 3 and 1.3 
mm respectively. The beam size of the telescope was 22 
for the J = 1 — » lines and 11 for the J = 2 — > 1 line. 
Pointing and focus were monitored regularly. The pointing 
corrections were never larger than 3 . The spectra were 
taken in position switching with a fixed reference position 
at (Z,&)=(0?65, 0?2), whic h was selected from the 13 CO 
map of Bally et al. ( 1987| ). Calibration of the data was 
made by observing hot and cold loads with known tem- 
peratures, and the line intensities were converted to main 
beam brightness temperatures, Tmb, using main beam ef- 
ficiencies of 0.68 and 0.41 for 3 and 1.3 mm respectively. 
The main beam efficiencies for the observations of June 
2000 are 0.80 and 0.53 for 3 and 1.3 mm respectively. 

A sample of spectra is shown in Fig. |2|. Most of the 
sources show CO emission in several velocity components 
with Gaussian line profiles. However, in some clouds the 
different components are blended, giving rise to more com- 
plex profiles. The observed parameters derived from Gaus- 
sian fits are listed in Table |^. 



2.2. ISO observations and results 

Several H 2 pure-rotational lines (from S(0) to S(5)) have 
also been observed towards the molecular clouds given in 
Table |l|. The observations were carried out with the Short 
Wavelength Spectrometer (SWS; de Graauw et al. 1996) 
on board ISO. The sizes of the SWS apertures at each 
wavelength are listed in Table ||. The orientation of the 
apertures on the sky varies from source to source, but 
it is within position angle 89.34° and 93.58° for all the 
observations (measuring the angles anti-clockwise between 
north and the short sides of the apertures). 

The observations presented in this paper are the re- 
sult of two different observing proposals. In one of them 
only the S(0), S(l) and S(3) lines were observed, in the 
second one all the lines from the S(0) to the S(5) but the 
S(2) were observed. The wavelength bands were scanned 
in the SWS02 mode with a typical on-target time of 
100 s. Three sources were also observed in the SWS01 
mode but the signal-to-noise ratio of these observations is 
rather poor and will not be discussed in this paper. Data 
were processed interactively at the MPE from the Stan- 
dard Processed Data (SPD) to the Auto Analysis Results 
(AAR) stage using calibration files of September 1997 and 
were reprocessed automatically through version 7.0 of the 
standard Off-Line Processing (OLP) routines to the AAR 
stage. The two reductions give similar results. In this pa- 
per we present the results of the reduction with OLP7.0. 
The analysis has been made using the ISAP2.O0 software 



The ISO Spectral Analysis Package (ISAP) is a joint de- 
velopment by the LWS and SWS Instrument Teams and Data 
Centers. Contributing institutes are CESR, IAS, IPAC, MPE, 
RAL and SRON. 
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package. With ISAP we have zapped the bad data points 
and averaged the two scan directions for each of the 12 
detectors. Then, we have shifted (flatfielded) the different 
detectors to a common level using the medium value as 
reference and finally, we have averaged the 12 detectors 
and rebinned to one fifth of the instrumental resolution. 
No defringing was necessary since the continuum flux at 
these wavelengths (A < 30/im) is lower than 30 Jy for all 
the clouds. 

Baseline (order 1) and Gaussian fitting to the lines 
have also been carried out with ISAP. The spectra are 
shown in Fig. || and the observed fluxes as derived from 
the fits are listed in Table ||. The absolute flux calibration 
errors are less than 30, 20, 25, 25, and 15% for the S(0), 
S(l), S(3), S(4), and S(5) lines, respectively (Salama et al. 



1997). Because of the medium spectral resolution of the 
SWS02 mode (A/AA - 1000 - 2000) and the wavelength 
calibration uncertainties (~ 15 — 50 km s -1 depending on 



the wavelength, see Valentijn et al. 1996 ), it is difficult to 
undertake a detail comparison between the kinematics of 
the H 2 lines and those of the 13 CO and C 18 lines. Table 
H lists the radial velocities of the S(l) lines, which have the 
higher signal-to-noise ratio. Within the calibration uncer- 
tainties, the radial velocity of the H2 lines agrees with at 
least one of the 13 CO components listed in Table 0. 

Unfortunately, the lack of resolution does not allow us 
to establish if the H2 emission is indeed arising from just 
one or several of the CO velocity components since, in 
general, all of them are within the velocity range of the 
unresolved H2 emission. M— 0.96 + 0.13 is the only cloud 
for which we can say that the warm H2 is not likely to arise 
in all the velocity components seen in CO. The CO com- 
ponents are centered at -110, 11, and 133 km s _1 , while 
the H2 S(l) line is centered at -70 km s _1 . Even with the 
spectral resolution of the SWS02 mode, one can see that 
the CO component with forbidden velocities (133 km s _1 ) 
is not likely to contribute to the H2 emission. 

Table|| also lists the widtns of the H 2 S(l) lines. The 

H2 line widths of the GC clouds tend to be larger than 

the instrumental resolution for extended sources (~ 170 
-1 



2000D . This 



km s" 1 for the S(l) line, see Lutz et al. 
is due to the large intrinsic line widths typical of the 
GC clouds and mainly, to the presence of several velocity 
components along the line of sight that contribute to the 
H2 emission. However, not all the sources that show CO 
emission in several velocity components have line widths 
larger than ~ 170 km s _1 (for instance M+0.83 — 0.10 or 
M+0.16 - 0.10). This implies that not all the CO veloc- 
ity components detected in these sources are contributing 
to the H2 emission, although it is difficult to discriminate 
which ones are emitting in H2. 

3. C ls O and 13 CO column densities 

The excitation analysis of the three lowest C 18 rota- 
tional lines by Hiittemeister et al. (1998) shows that the 



CO emission could arise in cold (20-30 K) and dense 
(10 4 cm' 3 ) gas or warmer (>100 K) and less dense gas 
(10 3 cm -3 ). However, the large column densities of cold 
(~ 25 K) dust (Martin-Pintado et al. 1999; Rodriguez- 
Fernandez et al. 2000) suggest that most of the J = 1 — » 
and J = 2 — > 1 CO emission should arise from cold and 
dense gas coupled to the dust. 

We have derived physical conditions and gas col- 
umn densities from the C 18 and 13 CO data using the 
Large Velocity Gradient (LVG) approximation (see e.g. 
Hiittemeister et al. 1998). Assuming a kinetic temperature 
(7k) of 20 K, we have constrained the H 2 densities (tih 2 ) 
from the J = 2^1toJ = 1^0 C 18 ratio (or the same 
ratio of 13 CO for a few sources, see Table |2|). Then, we have 
derived the C 18 and 13 CO column densities (N c is and 
Ni3qq) for corresponding nn 2 , using the J = 1 — > lines 
intensities. The results of the analysis are listed in Table 
|. Typical C 18 J = 2^1toJ=1^0 line ratios are ~ 
1.0-1.5 which give n Ha of - lO 3 - 5 " 4 cur 3 for T K =20 K. 
The typical integrated intensities of the C 18 0( J = 1 — >- 0) 



lines (~ 3-9 K km s ) imply Nq 



18Q 



of 



10 1 



cm . Ai3 CO is approximately a factor of 10 larger than 
Nqisq. Since the expected abundance ratio of the two 
species in the GC is 12.5 (Wilson & Matteucci [1992D , the 
observed ratio indicates that both lines are optically thin. 
We can also explain the observed CO lines ratios and in- 
tensities with higher kinetic temperatures (see above). For 
instance, for Tk=100 K one would obtain H2 densities 
which are lower by a factor of ~ 2.5. However, even in the 
unrealistic case that all the J = 1 — > and J = 2 — > 1 CO 
arise in warm gas with 7k=100 K, the column densities 
do not change more than 10% with respect to those at low 
temperature. Thus, in general, if one considers a mixture 
of warm and cool gas the total column densities traced by 
CO will be similar to those derived with 7k=20 K. 

Table || also gives the estimated H2 column densities 
as derived from Ai3co assuming that the abundance of 
13 CO relative to H2 is 5 10~ 6 . This ratio is based on the 
13 C/ 12 C isotopic ratio in the Galactic center of 1/20 (Wil- 
son & Matteucci |1992| ) and a CO/H 2 ratio of 10~ 4 (see 
e.g. Hiittemeister et al. 1998 and references therein). The 
typical H2 column densities derived for the main velocity 
components in all the sources are of a few 10 2 



4. Warm H 2 

Table || lists the observed fluxes of the H2 lines. The 
most intense lines are the S(0) ( J = 2 — > 0) and S(l) 
( J = 3 — * 1) lines, with typical fluxes of 0.5-1 10~ 19 and 
1-2 10~ 19 W cm -2 , respectively. Unfortunately, the S(2) 
line was only observed in the two clouds already discussed 
in detail by Rodriguez- Fernandez et al. ( pOOOp . The S(3) 
line is very weak and it has only been detected in the 
sources with more intense S(l) emission. Even in some 
sources which show emission in the S(4) and S(5) lines, 
the S(3) line has not been detected. This is due to strong 
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dust absorption produced by the solid state band of the 
silicates at 9.7/xm (Martin-Pintado et al. 199E). 

The pure rotational lines of H2 arise due to electric 
quadrupole transitions. The qua drupo le transition proba- 
bilities are small (Turner et al. 1977 ) and thus the rota- 
tional lines remain optically thin. In this case, the column 
density of the upper level involved in a transition from 
level i to level j can be obtained from the line fluxes Fij 
of Table || using the following expression: 



i 1 , , 

where A. 



i and 



Aij are the wavelength and the 
quadrupole probability of the transition, and Tij and 
are the dust opacity and the aperture at A^ , respectively. 
Since the column densities are averaged on the ISO aper- 
tures, in the case of extended sources (assumed homoge- 
neous) , it is not necessary to apply any additional correc- 
tion to account for the different ISO apertures (see also 
Rodriguez- Fernandez et al. 2000). 



4-1. Extinction and ortho-to-para ratio 

Figure ^ shows the population diagrams for one of the 
sources for which more than four lines were detected: 
M— 0.32 — 0.19. It shows, for each observed line, the log- 
arithm of the upper level population divided by both the 
rotational and nuclear spin degeneracy, i.e. 3(2 J + 1) for 
the ortho levels (odd J) and (2 J + 1) for the para levels 
(even J). 

The filled circles show the populations without any 
extinction correction. One can see that the population in 
the J=5 level is lower than expected from the interpo- 
lation from the other levels. As discussed in Rodriguez- 
Fernandez et al. (200C), this fact can be used to estimate 



the total extinction caused by the dust located between 
the observer and the H2 emitting region. Once an extinc- 
tion law is assumed, we can correct the H2 line intensities 
by increasing the visual extinction until the column den- 
sity in any level (in particular that in the J=5 level) is con- 
sistent with the column densities derived for the other lev- 
els, i.e. until the population diagrams are smooth curves. 

We have used the extinction law derived by 
Lutz ( |1999| ) towards the Galactic center using Hydrogen 
recombination lines. This extinction law differs from that 
of Draine ( |1989 ) for silicate-graphite mixtures of grains in 
that there is no deep minimum at ~ 7/im and there is a 
slightly higher value for the Ag^^/Ay ratio, where Ay 
the visual extinction (at 0.55 //m) and A).7 pm is the extinc- 
tion at 9.7 /Ltm. For instance, in the case of M— 0.32 — 0.19 
one sees that 15 mag of visual extinction (squares in Fig. 
||) is a lower limit to the extinction while 45 mag (stars in 
Fig. [|) is an upper limit. The best result is obtained for a 
visual extinction of around 30 mag (triangles). Using this 
method for the other sources with more than four lines 
detected, we also derive a visual extinction of ~ 30. This 
value should be considered as a lower limit to the actual 



extinction for the sources where the S(3) line was not de- 
tected. It is not not possible to know how much of this ex- 
tinction is caused by material in the line-of-sight towards 
the GC (foreground extinction) and how much is intrinsic 
to the GC clouds. Nevertheless, a visual extinction of ~ 30 
mag is in agreement with the average foreground extinc- 



tion as measured by Catchpole et al. (1990) using stars 
counts and suggests that the H2 emission can arise from 



the clouds surfaces (see also Pak, Jaffe & Keller 1996). In 



the other sources where we cannot estimate the extinction 
from our H2 data we have applied a correction of Ay=30 
mag. For those clouds located farther from the center of 
the Galaxy and/or the Galactic plane, we have corrected 
the observed fluxes by 15 mag (see Table |4|). This value 



was derived by Rodriguez-Fernandez et al. ( |2000| ) by an- 
alyzing the far infrared dust emission toward two sources 
in the "Clump 2" and the 1=1.5° complexes. In any case, 
the extinction correction has a minor impact in the main 
results of this paper (see below). Figure || shows the ex- 
tinction corrected population diagrams for all the sources 
presented in this paper. 

The values of extinction required to give a smooth 
population diagram would be somewhat smaller if the H2 
ortho-to-para (OTP) ratio were lower than the local ther- 
modynamic equilibrium (LTE) ratio. This is obvious since 
the method to derive the extinction depends mainly on 
the extinction at the wavelength of an ortho level (J=5). 
Non-equilibrium OTP ratios measured with the lowest ro- 
tational lines has been found in two clouds of our sample 
(Rodiguez-Fernandez et al. 2000). Unfortunately, for the 
clouds presented in this paper, it is difficult to estimate the 
OTP ratio since the S(2) line has not been observed and 
the S(3) line is completely extincted in most of them. Cur- 
rent data do not show any evidence for a non-equilibrium 
OTP ratio, but we cannot rule it out a priori. For in- 
stance, assuming OTP ratios of ~ 2 we still can find a 
smooth population diagrams, i.e. without the typical zig- 
zag shape characteristic of non-equilibrium OTP ratios 
(see. e.g. Fuente et al. 1999). In this case, the extinction 
would be of ~ 20 — 25 mag instead of 30 mag. On the 
contrary, assuming OTP ratios of ~ 1 one finds, in gen- 
eral, rather artificial diagrams, which suggests that OTP 
ratios as low as ~ 1 are not compatible with the data. 
Although one must bear in mind these considerations, in 
the following we assume that the OTP ratios are LTE. 

4-2. H2 column densities and excitation temperatures 

Table [j] lists the results derived from the H2 lines after ap- 
plying the extinction corrections. The excitation tempera- 
ture derived from the J=3 and J=2 levels (T 32 ) is between 
130 and 200 K while the excitation temperature derived 
from the J=7 and J=6 levels (T 76 ) is - 500 - 700 K. 
The temperatures are only 15-20 % larger than those one 
obtains without any extinction correction. For the four 
clouds in which the S(4) and S(5) lines were undetected, 
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we derive T 32 of ~ 135-150 K, clearly lower than the tem- 
perature derived for the sources where the S(4) and S(5) 
lines were detected. There is no clear dependence of T 32 on 
the distance to the Galactic center. However, it is notewor- 
thy that the two clouds with lower T 32 are located in the 
Sgr C complex, one of them in the non-thermal filament. 

Obviously, T 32 lacks of physical sense if the ortho-H 2 
and para-H 2 abundances are not in equilibrium. As men- 
tioned, we can obtain smooth population diagrams assum- 
ing OTP ratios lower than the LTE ratio. The tempera- 
ture T 32 derived in this case (T^™) is higher than the one 
derived directly from the observations (T 32 ). For instance, 
assuming OTP ratios ~ 2 one obtains a T 32 rr which is 
- 10% larger than T 32 . 

It is possible to estimate the total warm H 2 column 
densities (Nn 2 ) by extrapolating the populations in the 
J=2 level to the J=l and J=0 levels at the tempera- 
ture T 32 . The derived warm 7Vh 2 are listed in Table |I] and 
should be considered lower limits to the actual amount of 
warm molecular gas since the lowest levels can be pop- 
ulated with colder, although still warm, gas. The total 
column density of warm H 2 varies from source to source 
but it is typically of 1-2 10 22 cm -2 . These column den- 
sities are only a factor of 1.2 higher than those one would 
obtain without any extinction correction. Thus, in regard 
to the derived gas temperatures and total column densi- 
ties, the extinction correction is not critical. On the other 
hand, extrapolating the column densities in the J=6 and 
J =7 to lower levels at the temperature T76, one finds that 
the amount of gas at ~ 600 K is less than 1% of the col- 
umn densities measured at ~ 150 K. The H 2 total column 
densities at temperatures T , | 2 rr assuming an OTP ratio of 
~ 2 are lower than those of Table | by a factor of 1.8. 
Note, that in this case the total column density should be 
derived extrapolating the observed population in the J=3 
to the J=l level and the population in the J=2 to the 
J=0 levels, as two different species at temperature T 32 rr . 
Of course, these column densities are still lower limits to 
the actual warm H 2 column densities. 

These results are the first direct estimation of the H 2 
column densities and the structure of the warm gas in 
the GC clouds. They show the presence of large column 
densities of warm molecular gas with large temperature 
gradients (150-700 K), ex tending the results derived by 



Huttemeister et al. ( 1993 ) from their NH 3 data. 



5. Discussion 

5.1. Warm H 2 to CO and NH3 ratios 

As mentioned in Sect. 2.2, we cannot identify which veloc- 
ity components seen in CO are associated to the warm H 2 . 
Furthermore, the bulk of the CO seen in the J — 1 — > 
and J — 2 — > 1 lines do not show the characteristics of 
warm CO associated to the warm H 2 (see Sect. 3). In the 
following, we will estimate the ratio of the warm H 2 col- 



umn densities observed with ISO to the H 2 column densi- 
ties derived from the CO using LVG calculations. We have 
added the column densities of each velocity component in 
every source. These total H 2 column densities are listed in 
Table § One can compare the maximum column densities 
derived from CO to the H 2 column densities listed in Table 
^ to derive a lower limit to the fraction of warm molec- 
ular gas with respect to the gas emitting in CO. These 
ratios are given in Table || for all the molecular clouds. 
We find that the warm H 2 is about ~ 30% of the H 2 col- 
umn densities measured from CO. For a few clouds the 
fraction of warm gas is as high as 77% (M— 0.50 — 0.03) 
or even ~ 100% for M-0.96 + 0.13. This implies that, for 
two clouds all the CO emission should arise from warm 
gas. 

Table || also lists the NH 3 abundances in the warm 
(X(NH 3 ) warm ) and cold components (X(NH 3 ) co id). The 
A(NH 3 ) warm has been derived from the column densi- 
ties of warm ammonia (Huttemeister et al. 1993) and our 
warm H 2 column densities. We find that, A(NH 3 ) warm is 
within a range of 3 10~ 8 to 4 10" 7 . On the other hand, 
A(NH 3 ) co id has been derived from the cold ammonia col- 
umn densities of Huttemeister et al. and the H 2 column 
densities derived from the 13 CO data. In this case, we have 
taken into account only the 13 CO velocity components 
with NH 3 emission and we have assumed that, in average, 
~ 70% of the gas traced by CO is cold gas. With these as- 
sumptions, ^(NH 3 ) co id varies between 4 10~ 8 and 6 10~ 6 , 
being the average value ~ 5 10~ 7 . This is similar to the 
abundance in the warm component, and approximately 
10 times higher than t he "ty pical" interstellar ammonia 
abundance (Irvine et al. 1987). The high NH 3 abundances 
in the cold gas point to the existence of a cold post-shocked 
gas component as suggested by Huttemeister et al. (1998) 
to explain the SiO emission in the GC clouds. 

5.2. Heating mechanism 

What is the heating mechanism that produces such a 
large amount of warm molecular gas in the GC? Shocks 
have been invoked to explain the widespread distribution 
and the large abundances of refractory molecules like SiO 
(Martin-Pintado et al. 1997, Huttemeister et al. 1998), the 
high temperatures observed in NH 3 (Wilson et al. 1982, 
Giisten et al. 1985) and the non-equilibrium H 2 ortho- 
to-para ratio of two sources in our sample (Rodriguez- 
Fernandez et al. 2000). The high NH 3 abundance derived 
in the previous section points to a mechanical heating 
mechanism since the ammonia molecule is easily photo- 
dissociated by ultraviolet radiation. The small column 
densities of warm dust in these clouds also points to a me- 
chanical heating mechanism (Martin-Pintado et al. 1999). 

On the other hand, in some of the clouds we have de- 
tected line emission from ionized species like Ne II, Ne III or 
O in, that should arise in an H 11 region ionized by ultravio- 
let (UV) photons (Martin-Pintado et al. 1999, 2000). This 
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implies that, at least in those clouds, there must be a PDR 
in the interface between the H 11 region and the molecular 
material. Large scale emission of the H2 u=l-0 S(l) line 
has also been interpreted as arising from PDRs of density 
n ~ 10 4 cm" 3 and incident far-UV flux of G ~ 10 3 (in 
units of 1.6 10" 3 ergs cm -2 s" 1 ) in the clouds surfaces 
(Pak, Jaffe & Keller 1996). The total visual extinction of 
~ 30 mag derived for the clouds of our sample matches 
the expected foreground extinction and suggest that the 
pure-rotational H2 emission could also arise in the surfaces 
of the clouds as the ro- vibrational lines. 

We have compared the population diagrams obtained 
for the GC clouds with the same type of diagrams pre- 
dicted by models of C-shocks, J-Shocks and PDRs. Fig- 
ure ^a shows the comparison between the predictions of 
a C-Shock from Draine et al. (1983), a J-Shock from Hol- 
lenbach & McKee ( |1989| ), and the data for M-0.32-0.19. 
The S(l) to S(5) lines (squares) can be explained with 
both a C-Shock with velocity of ~ 12 km s" 1 acting on 
gas with preshock density of 10 6 cm" 3 (circles) or a J- 
Shock of 50 km s" 1 and preshock density of 10 6 cm -3 
(triangles). However, the observed emission in the S(0) line 
is ~ 3 times larger than the predicted by both models. 

Figure ^b shows the population diagram for M+0.16 — 
0.10 (squares) versus the prototypical reflection nebula 
NGC 7023 (triangles). As discussed by Fuente et al. 
(1999), the H2 emission from this so urce i s well fitted 
by the PDR model of Burton et al. ( |1990| , |1992[ ) with 



G = 10 4 and n=10 6 cm" 3 although with an OTP ratio 
of 1.5-2. Comparing the NGC 7023 population diagram 
with M+0.16 — 0.10, one finds that the agreement is ex- 
cellent for the S(4) and S(5) lines but it is not so good 
for the lowest lines, even taking into account the non- 
equilibrium OTP ratio found in NGC 7023. In particular, 
the GC clouds exhibit more emission in the lowest lines 
than expected from the PDR model for Go = 10 4 and 
n=l0 6 cm" 3 . In contrast, the H2 v=l-0 S(l) intensity 
predicted by this PDR model is a factor of ~ 10 larger 
than observed by Pak, Jaffe & Keller (1996). This fact 
would imply that the vibrational line emission is more di- 
luted in a 3 beam than the pure-rotational lines in the 
SWS beam or that the PDR models do not apply. 

In any case, the observed curvature of the population 
diagrams seems to be in good agreement with the pre- 
dicted temperature gradient in a PDR. In Fig. |d, we also 
show the population diagram one obtains integrating the 
H2 emission in LTE with the temperature and H2 abun- 
dance profiles along the Go=10 4 and n=10 6 cm" 3 PDR 
model of Burton et al. (1990). The result differs from that 
of Burton et al. in that we do not take into account any 
radiative pumping, which affects mainly to higher levels 
than those involved in the S(0) and S(l) lines. Althoght 
the GC emission is ~ 3 times larger, it is evident that the 
shape of the population diagram is very similar to that 
observed. 



With regard to those sources where the S(4) and S(5) 
were not detected, the upper limits imply that if they 
are PDR-excited the density must be somewhat lower 
than rt=10 6 cm" 3 , or if shock-excited, the shock velocity 
should be slightly lower than those of the models ploted 
in Fig. |. 

Both shock and PDR models suggest densities as high 
as 10 6 cm -3 and fail to explain the observed intensity of 
the S(0) emission and to less extend the S(l) line. The den- 
sities implied by the models seem somewhat large, but it 
looks like the H2 traces two components: a hot (~ 500 K) 
and dense (<10 6 cm" 3 ) component necessary to explain 
the observed S(4) and S(5) lines, and a warm component 
(~ 150 K) traced by the S(0) and S(l) lines. To match 
the measured J = 2-^l/J=1^0 13 CO and C 18 ra- 
tios the warm H2 component should have densities of 
~ 10 3 cm" 3 (see Sect. 3). The hot and dense gas would 
have J = 2— >1/J=1— >0 13 CO ratios of - 4 - 5 but it 
would emit mainly in the high-J CO lines. In any case, 
the column density of hot and dense gas is very small to 
make it detectable in the low- J CO lines when mixed with 
the colder and less dense gas that dominates the emission 
of these lines. 

To explain the derived T32 ~ 150 K is necessary to 



invoke PDRs with Go ~ 10 and n 



10 3 



but 



to obtain the observed intensities ~ 20 of such PDRs are 
needed. J-shock models do not predict temperatures as 
low as 150 K. Moreover, the high velocities required to 
explain our data are difficult to reconcile with the obser- 
vations. C-shocks could explain the observed S(0) and S(l) 
emission with, at least, 10 shocks with velocities as low as 
~ 7 km s" 1 and n=10 6 cm" 3 (even more shock fronts are 
needed for lower gas densities). In addition, dissipation of 
supersonic turbulence could heat the gas to temperatures 
of - 150 K (Wilson et al. |1982| ; Giisten et al. p.985| ) and 
thus, could contribute to the emission in the two lowest H2 
lines. The origin of the turbulence would be the movement 
of dense clumps in a less dense inter-clump medium due to 
the differential Galactic rotation and the tidal disruption 
of the clumps. 

The heating rate by dissipation of supersonic turbu- 
lence can be estimated as T ~ 3.5 10 28 v 3 nH 2 (lpc/Z) erg 



s 1 cm 3 (Black 1987| ), where I and vt are the spa- 



tial scale and the velocity of the turbulence, respectively. 
Taking vt ~ 15 km s" 1 (the typical linewidths of GC 
clouds), I — 5 pc 

r 

gas, T ~ 150 K and nn 2 ~ 10 3 cm" 3 , the cooling is 
expected to be dominated by H2 and CO. Le Bourlot 



and tih 2 = 10 3 cm 3 , one obtains 
5 10~ 22 erg s" 1 cm" 3 . For the conditions of the warm 



et al. (1999) has recently estimated the cooling rate by 
H2 (Ah 2 ) lor a wide range of parameters. For the warm 
gas component of the GC clouds we obtain Ah 2 ~ 3 10 -22 
erg s _1 cm -3 , which is comp arable to the CO cooling rate 



(see e.g. Goldsmith & Langer 1978 ) Thus, comparing heat- 
ing and cooling rates, one finds that the dissipation of su- 
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personic turbulence could account for the heating of the 
warm component. 

In summary, several agents could heat the warm com- 
ponent, while the hot component should trace the densest 
gas in the GC clouds heated by a PDR or a shock. For in- 
stance, if the inhomogeneous structure revealed in the Sgr 
B2 envelope by interferometric NH3 observations (Martin- 
Pintado et al. 1999) is common in the GC, and due to 
evolved massive stars as they propose, both C-shocks of 
~ 10 km s _1 (shell expansion) and PDRs (stellar radia- 
tion) would be present. However, it is not possible to rule 
out mechanical heating by large scale shocks. In fact, the 
high fraction of warm H2 derived for M— 0.96 + 0.13 and 
the fact that the CO component with positive velocities 
apparently does not contribute to the H 2 emission sug- 
gests this kind of heating since, at this galactic longitude, 
shocks are expected at negative velocities due to the in- 
tersection of X\ and xi orbits in the context of a barred 



potential (Binney et al. 1991). 



6. Summary and conclusions 

We have observed the S(0) to S(5) H 2 pure-rotational lines 
with the SWS spectrometer on-board ISO toward a sam- 
ple of 18 molecular clouds of the Galactic center region. 
The S(3) line is strongly affected by dust extinction due 
to the 9.7 /im band of the silicates. After correcting the 
H2 data for extinction using a self-consistent method, and 
assuming that the ortho- and para-Bi2 populations are in 
equilibrium one finds that the S(0) and S(l) lines indi- 
cate temperatures of ~ 150 K. Extrapolating to the lowest 
levels at that temperature, a total H2 column density of 
~ 1 — 2 10 22 cm -2 is derived. This is the first direct esti- 
mate of the column density of warm gas in the GC clouds. 
In addition, it shows a complex temperature structure of 
the warm gas. 

The temperature derived from the S(5) and S(4) levels 
is ~ 600 K in the sources in which it can be derived, how- 
ever the column density of gas at this temperature is less 
than 1% of the column density at T=150 K. Assuming an 
OTP ratio of ~ 2 the temperatures would be 10% larger 
than those derived assuming a LTE OTP ratio, while the 
total H2 column densities at those temperatures would be 
a factor of ~ 1.8 lower than the column densities derived 
assuming the ortho- and para-Bi2 populations in equilib- 
rium. Comparing the H2 warm column densities with the 
column densities derived from our CO data by LVG calcu- 
lations one finds that the average fraction of warm H2 to 
the gas observed in CO is ~ 30%. With our data and the 
NH3 observations of Huttemeister et al. (1993) we derive 
relatively high NH3 abundances of a few 10~ 7 in both the 
warm and the cold components. 

Several indirect arguments point to shocks as the heat- 
ing mechanism of the warm gas but PDRs may also play 
a role. Direct comparison of the H2 data with PDRs and 
shocks models indicate that the S(4) and S(5) trace the 



densest gas in the GC clouds (<10 6 cm~ 3 ) heated in PDRs 
or shocks. Nevertheless, such dense PDRs or shocks fail to 
explain the S(0) and S(l) lines: several less dense PDRs, 
low velocity shocks (< 10 km s" 1 ) or both along the line 
of sight would be needed to explain the observed emission. 

The cooling by H2 in the warm component of GC 
clouds is comparable to the cooling by CO. Equating the 
H2 cooling rate with the heating rate by dissipation of su- 
personic turbulence, one finds that this mechanism could 
also contribute to the emission in the two lowest H2 lines. 
In one source (M— 0.96 + 0.13), we have also found some 
evidence of large scale shocks that should be checked with 
higher spectral resolution H2 observations. 
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Table 1. J2000 coordinates of the sources 



Source 


RA 


DEC 


Complex 




h m s 






M-0.96 + 0.13 


17:42:48.3 


-29:41:09.1 


SgrE 


M-0.55 - 0.05 


17:44:31.3 


-29:25:44.6 


SgrC 


M-0.50 -0.03 


17:44:32.4 


-29:22:41.5 


SgrC 


M-0.42 + 0.01 


17:44:35.2 


-29:17:05.4 


SgrC 


M-0.32 -0.19 


17:45:35.8 


-29:18:29.9 


SgrC 


M-0.15 -0.07 


17:45:32.0 


-29:06:02.2 


Sgr A 


M+0.16 -0.10 


17:46:24.9 


-28:51:00.0 


Arc 


M+0.21 -0.12 


17:46:34.9 


-28:49:00.0 


Arc 


M+0.24 + 0.02 


17:46:07.9 


-28:43:21.5 


Dust Ridge 


M+0.35 - 0.06 


17:46:40.0 


-28:40:00.0 




M+0.48 + 0.03 


17:46:39.9 


-28:30:29.2 


Dust Ridge 


M+0.58 -0.13 


17:47:29.9 


-28:30:30.0 


Sgr B 


M+0.76 - 0.05 


17:47:36.8 


-28:18:31.1 


SgrB 


M+0.83-0.10 


17:47:57.9 


-28:16:48.5 


SgrB 


M+0.94 - 0.36 


17:49:13.2 


-28:19:13.0 


SgrD 


M+2.99 - 0.06 


17:52:47.6 


-26:24:25.3 


Clump 2 





1 1 1 1 | 1 1 1 


1 | ■ 1 1 1 | I 1 I 1 | 1 I 1 l_ 


45 




M + 0.32-0.19 — 




~ J = 3\V 




40 












35 
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E(K) 

Fig. 4. Population diagrams for M— 0.32 — 0.19 without 
any extinction correction (circles) and corrected for 15 
(squares), 30 (triangles), and 45 mag. (stars) of visual ex- 
tinction {A\). We have assumed the relative extinctions 
derived toward the Galactic center by Lutz (1999). Note 
that a smooth curve is obtained with Ay ~ 30 mag. 
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Table 4. Total H2 column densities and rotational tem- 
peratures between the J=3 and the J =2 levels (T32) and 
between the J=7 and the J=6 levels (Tre) after correcting 
for extinction. Numbers in parentheses are la errors of the 
last significant digits as derived from the fluxes errors in 
the Gaussian fits of the lines. 



Source 


Ay 


T32 


Tre 








K 


K 


10 22 cm -2 


M-0.96 + 0.13 


15 


157(6) 




1.10(9) 


M-0.55 -0.05 


30 


135(5) 




2.7(3) 


M-0.50 -0.03 


30 


135(4) 




2.3(2) 


M-0.42 + 0.01 


30 


167(6) 




1.03(8) 


M-0.32 - 0.19 


30 


188(5) 


650(90) 


1.03(5) 


M-0.15 - 0.07 


30 


136(6) 




2.6(4) 


M+0.16 - 0.10 


30 


157(7) 


900(200) 


1.17(13) 


M+0.21 -0.12 


30 


186(13) 


670(110) 


0.64(7) 


M+0.24 + 0.02 


30 


163(2) 




1.73(6) 


M+0.35 - 0.06 


30 


195(11) 


700(200) 


0.66(5) 


M+0.48 + 0.03 


30 


174(7) 


<600 


1.03(9) 


M+0.58 -0.13 


30 


149(5) 




1.3(2) 


M+0.76 - 0.05 


30 


181(4) 




1.77(8) 


M+0.83 -0.10 


30 


178(5) 


550(60) 


1.59(6) 


M+0.94 - 0.36 


15 


146(7) 




0.95(10) 


M+2.99 - 0.06 


15 


152(3) 




1.40(9) 



45 



35 



45 



I 


1 1 1 1 1 1 1 1 1 1 1 1 
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Fig. 6. a: Population diagram for M— 0.32 — 0.19 (open 
squares) corrected for 30 mag. of visual extinction. The 
errorbars represent upper limits to the flux calibration un- 
certainties (see text). For comparison, it also displays the 
population diagrams derived from the model of Draine 
et al. (1983) of a shock with velocity ~ 12 km s _1 and 
preshock density 10 6 cm^ 3 (circles and dashed lines). Tri- 
angles and long-dashed lines are used to plot the popula- 
tion diagram derived from the J-shock model of Hollcn- 
bach & McKee (1989) for a velocity of 50 km s _1 and a 
preshock density of 10 6 cm~ 3 . b: Comparison of the pop- 
ulation diagram derived for M+0.16 — 0.10 (open squares) 
with the results of Fuente et al. (1999) for the NGC 7023 
PDR (triangles and dashed line) and the population di- 
agram one obtains integrating the H2 emission along the 
temperature and H2 abundance gradient derived by Bur- 
ton et al. (1990) for a PDR with density of 10 6 cm" 3 and 
Go = 10 4 (open circles) 
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Table 2. Observational parameters and LVG results for the CO data: Integrated intensities of the J = 1 — > transitions 
of C 18 and 13 CO and C 18 J = 2 ^ 1 to J = 1 ^ line intensity ratio. Column densities and nn 2 derived from 
the LVG calculations. 7Vh 2 derived from Niz C q assuming a 13 CO abundance relative to H2 of 5 10~ 6 . Numbers in 
parentheses are la errors of the last significant digit. 



Source 


13 co 
km s -1 


r 13 co 
(1-0) 
K km s -1 


rC ls O 

K km s _1 


tC 1s O 1 rC 18 
- , (2-l)/ J (l-0) 


'op(n H2 ) a 
log( cm -3 ) 


7V cl80 a 
10 15 cm" 2 


iVl3 CO a 

10 16 cm" 2 


Nh 2 
10 21 cm" 2 


M-0.96 + 0.13 


-110 


9.9(7) 




1.2(3) C 


3.5-4 




0.54-1.1 


1.1-2.2 




11 


24(1) 




0.7(2) c 


3.2-3.6 




0.54-2.1 


1.1-4.2 




133 


30.2(6) 




0.8(l) c 


3.3-3.5 




1.9-2.4 


3.8-4.8 


M-0.55-0.05 


-102 


126(3) 


11.2(5) 


2.2(1) 


4-4.4 


8.9-14.1 


14.4-18.2 


28.8-36.4 




-60 


45(2) 


3.2(4) 


2.5(4) 


>4.2 


3.1-5.0 


5.4-8.6 


10.8-17.2 




-2 


23(2) 


3.5(4) 


1.7(3) 


3.8-4 


2.7-3.3 


1.9-3.0 


3.8-6.0 


M-0.50-0.03 


-98 


119(5) 


20(2) 


0.9(1) 


3.4-3.7 


12-15 


12-15 


24-30 


M-0.42 + 0.01 


-87 


66(4) 


6.6(11) 


1.7(4) 


4-4.5 


4.5-7.1 


6.0-9.5 


12-19 




-74 


54(4) 


3.0(8) 


2.0(7) 


>4 


1.9-3.8 


4.7-7.5 


9.4-15 


M-0.32-0.19 b 


-77 


16(2) 






>3 




0.95-2.4 


1.9-4.8 




-31 


40(2) 






>3 




2.7-5.4 


6.4-10.8 




23 


24(2) 






>3 




1.5-3.4 


3.0-6.8 


M-0.15-0.07 


5 


120(7) 


16.1(4) 


1.3(1) 


3.7-3.9 


1.3-1.6 


20-25 


40-50 




36 


155(7) 


13.8(7) 


1.2(1) 


3.5-4.1 


8.3-13.1 


13-17 


26-34 


M+0.16 - 0.10 


-8 


38(5) 


2.8(4) 


1.7(3) 


3.8-4.2 


1.8-2.8 


3.2-4.5 


6.4-9 




04 


97(8) 


9.2(5) 


2.0(2) 


3.9-4.1 


6.9-8.7 


12.7-16 


25.4-32 




81 


34(8) 


2.7(4) 


1.6(4) 


>3.7 


1.6-3.5 


2.5-4.0 


5-8 


M+0.21 - 0.12 b 


66 


42(3) 




0.3(1) 


>3 




2.9-5.9 


5.8-11.8 




100 


14(3) 




0.2(1) 


>3 




0.92-1.8 


1.8-3.6 


M+0.24 + 0.02 


30 


190(5) 


17.7(11) 


1.0(1) 


3.5-3.9 


10-16 


16-25 


32-50 




66 


66(3) 


7.0(5) 


1.1(2) 


3.6-3.8 


4.5-5.6 


8.2-10.3 


16.4-20.6 


M+0.35 - 0.06 


22 


44(2) 




1.7(l) c 


4-4.1 




3.8-6.0 


7.6-12.0 




75 


38.7(13) 




1.9(l) c 


4.4-4.8 




3.1-4.9 


6.2-9.8 




94 


24.2(11) 




1.4(1)° 


3.8-4 




1.7-2.5 


3.4-5.0 


M+0.48 + 0.03 


14 


58(7) 




1.4(2) c 


3.8-4.2 




4.8-6.0 


9.6-12.0 




33 


107(7) 




l.l(l) c 


3.6-3.8 




11-12 


22-24 


M+0.58-0.13 


36 


143(4) 




l.l(l) c 


3.6-3.7 




13-16.7 


26-33.4 




87 


36(4) 




1.3(2) c 


3.7-3.9 




2.4-3.0 


4.8-6.0 


M+0.76 - 0.05 


-40 


20(3) 


3.4(11) 


0.6(4) 


2.8-4.2 


1.8-4.4 


1.3-5.2 


2.6-10 




33 


281(5) 


31(2) 


1.6(2) 


3.7-4.3 


19-37 


25-30 


50-60 




84 


72(3) 


8.0(11) 


0.9(2) 


3.3-3.6 


4.1-6.5 


6.5-8.2 


13-16 


M+0.83 - 0.10 


-40 


9.4(9) 


1.8(4) 


<0.8 


<3.2 


>0.9 


>5.6 


>11.2 




5 


46(3) 


6.1(7) 


0.6(1) 


3.0-3.3 


3.4-6.8 


3.4-5.4 


6.8-10.8 




32 


93(3) 


7.7(7) 


1.4(2) 


3.6-4.5 


4.5-9.2 


7.5-12 


15-24 




96 


88(2) 


9.2(5) 


1.3(1) 


3.6-3.8 


5.6-7.1 


7.6-9.6 


15.2-19.2 


M+0.94 - 0.36 b 


-83 


13.8(11) 






>3 




0.9-1.8 


1.8-3.6 




-67 


18.1(12) 






>3 




1.2-2.3 


2.4-7.2 




-36 


39(2) 






>3 




2.6-5.2 


5.2-10.4 




-10 


28.2(14) 






>3 




1.9-3.7 


3.8-7.4 


M+2.99 - 0.06 


-12 


37(4) 


2.6(14) 


2.1(14) 


>3.7 


1.5-3.0 


2.4-4.7 


4.8-9.4 




12 


34(4) 


3.5(11) 


0.7(4) 


3.0-4.0 


2.4-3.8 


2.4-6.0 


4.8-12.0 



a Minimum and maximum values derived with the LVG code taking into account the la errors in the line intensi ties 
b When the ratio 2-1/1-0 could not be obtained, we have assumed riH 2 > 10 3 cm" 3 following Hiittemeister et al 



CW-o° instead of I^/I^ 



1998. 
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Table 3. Fluxes of the H2 lines as derived from Gaussian fits in units of 10 -20 W cm -2 . Upper limits are 3a values 
at the instrumental spectral resolution for point sources. Numbers in parentheses are la errors of the last significant 
digit as derived from the Gaussian fits. The radial velocities and the widths of the lines with better signal-to-noise 
ratio (the S(l) lines) are also shown. The errors in the radial velocities are dominated by the wavelength calibration 
uncertainties (15-30 km s^ 1 for the S(l) line). Typical la error of the line widths derived from the Gaussian fits is 
less than 5 km s^ 1 . 



Line 




") 


S(0) 


S(l) 


S(3) 


S(4) 


S(5) 






Aper. ( 


X 


20 x 27 


14 x 27 


14 x 20 


14 x 20 


14 x 20 


km s _1 


km s 


A(^/m) 






28.2188 


17.03483 


9.66491 


8.02505 


6.9095 






M-0.96 


+ 


13 


7.8(9) 


18.4(8) 


2.2(5) 






-70 


270 


M-0.55 


-0 


05 


9.5(14) 


9.7(6) 


<0.80 


<0.78 


<2.0 


-80 


230 


M-0.50 


- 


03 


8.2(10) 


8.4(4) 


<0.64 






-60 


230 


M-0.42 


+ 


01 


6.2(6) 


13.1(7) 


<0.70 






-57 


230 


M-0.32 


- 


19 


7.8(6) 


23.0(6) 


2.1(2) 


3.5(7) 


5.7(8) 


-59 


230 


M-0.15 


- 


07 


9.4(13) 


9.9(12) 


<1.1 


<1.4 


<2.8 


-35 


220 


M+0.16 


-0 


10 


6.1(9) 


10.5(7) 


<0.9 


2.7(6) 


6.5(10) 


40 


180 


M+0.21 


-0 


12 


4.7(9) 


13.3(8) 


<1.2 


2.8(4) a 


4.8(ll) a 


16 


260 


M+0.24 


+ 


02 


9.8(5) 


18.9(4) 


<0.92 






-6 


170 


M+0.35 


-0 


06 


5.3(8) 


17.2(6) 


<1.0 


2.0(7) 


3.5 (8) 


27 


200 


M+0.48 


+ 


03 


6.7(8) 


15.9(8) 


1.6(3) a 


2.4(10) a 


<3.4 


17 


170 


M+0.58 


-0 


13 


6.0(6) 


8.7(7) 


<0.96 


<0.97 


<2.1 


4 


210 


M+0.76 


- 


05 


12.4(9) 


32.8(8) 


2.0(5) 






-18 


180 


M+0.83 


-0 


10 


10.8(9) 


27.1(4) 


2.2(3) 


5.6(10) 


6.7(8) 


16 


170 


M+0.94 


- 


36 


5.7(9) 


10.6(5) 


<1.2 


<1.1 


<2.7 


-30 


190 


M+2.99 


-0 


06 


9.2(6) 


19.4(8) 


<0.79 






28 


190 



a Detections with low signal-to-noise ratio (~ 2.5) 




Fig. 1. The positions of all the sources of our sample (including the two clouds presented in Rodriguez- Fernandez 
et al. |2000| ) overlayed in the C 18 O(l-0) map by Dahmen et al. (1997). 
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M-0.96+0.13 M-0.55-0.05 M+0.35-0.06 M+0.83-0.10 




Fig. 2. 13 CO and C 18 spectra of four sources. 

Table 5. Total column densities of H2 derived from 13 CO. Fraction of warm Nn 2 as measured with ISO to the total 
-/Vh 2 derived from 13 CO. Abundances of NH 3 in the warm and cold components (see text). 



Source 


Nu 2 C0 
10 22 cm -2 


jV H2 warm /iVH 2 uu 


X(NH3) warm 


X(NH 3 ) co id 


M-0.96 + 0.13 


0.6-1.1 


1 


3.7 10- y 


4.9 10~ b 


M-0.55 - 0.05 


4.3-6.0 


0.45 






M-0.50 - 0.03 


2.4-3.0 


0.77 


2.6 10~ 8 


1.6 10~ 7 


M-0.42 + 0.01 


2.1-3.4 


0.29 


8.3 10~ 8 


2.9 10~ 8 


M-0.32 - 0.19 


1.1-2.2 


0.45 


1.8 10~ 8 


3.1 10~ 7 


M-0.15-0.07 


6.6-8.4 


0.31 


2.4 10~ 7 


2.7 10~ 7 


M+0.16 - 0.10 


3.7-4.9 


0.24 






M+0.21 - 0.12 


0.8-1.5 


0.41 






M+0.24 + 0.02 


4.8-7.1 


0.24 


1.3 10~ 7 


8.9 10~ 7 


M+0.35 - 0.06 


1.7-2.7 


0.25 






M+0.48 + 0.03 


3.2-3.6 


0.28 






M+0.58-0.13 


3.1-3.9 


0.33 






M+0.76 - 0.05 


6.6-8.6 


0.21 






M+0.83-0.10 


4.8-6.5 


0.25 




3.4 10~ 8 


M+0.94 - 0.36 


1.3-2.9 


0.33 




6.7 10~ 7 


M+2.99 - 0.06 


1.0-2.1 


0.65 




9.0 10~ 7 
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Fig. 5. Population diagrams for all the sources corrected for the extinctions listed in Table |4| The filled circles are 
connected when more than three lines are detected. Arrows indicate upper limits. The error-bars are smaller than the 
circles (even taking into account both calibration and Gaussian fitting errors). 



